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Abstract: The low-energy excited states of a system of interacting one-dimensional 



\0 I fermions in a conducting state are collective charge and spin density oscillations. 

The unusual physical properties of such a system (called "Luttinger liquid") are 
CN ■ characterized by the velocities Up and Mo- of the charge and spin excitations, as well as 

by a parameter Kp that determines the power law behavior of correlation functions. 
Q"^ I Umklapp scattering occuring at half-filling or other commensurate bandfilling can 

lead to a transition into an insulating state, characterized in particular by a gap in 
the charge excitations (the Mott-Hubbard gap). The properties in the vicinity of 
the transition are shown to depend on both the way the transition is approached 
(constant bandfilling and varying interaction, or constant interaction and varying 



Q I bandfilling) and on the "order" of the commensurability. In particular, even and odd 

fractional fillings show quite different behavior. This behavior is illustrated in detail 
using lattice models like the Hubbard model and its extensions. 
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■ 1 Introduction 

A theoretical understanding of interacting fermion systems in one dimension is impor- 
tant for a number of reasons. On the one hand, in the physics of quasi-one-dimensional 
organic conductors or of conducting polymers interaction effects play a major 
role. On the other hand, one-dimensional models can be easier to understand than 
their higher- dimensional versions, or even exactly solvable, as is the case with the 
prototypical model of correlated fermions, the Hubbard model 0. They therefore 
can provide valuable information on the role of correlation effects in higher dimension, 
e.g. on the physics of correlated fermions in two dimensions which is thought to be 
at the origin of the many interesting properties of high-temperature superconductors 

Theoretical work on interacting fermions in one dimension has progressed along 
a number of different lines. One approach has been the perturbative investigation 
of the weak coupling limit. Even this is in fact not entirely straightforward, mainly 
because of the infrared divergences encountered in this type of calculation which re- 
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quire a renormalization group treatment. A complete review of this approach has 
been given by Solyom An alternative and more general approach is provided by 
the so-called "bosonization" method, which is based on the equivalence between in- 
teracting fermions and noninteracting bosons (representing density fluctuations) and 
on the expression of fermionic operators in terms of these bosons. Combined with the 
renormalization group approach, the bosonization method provides a rather straight- 
forward description of the peculiar properties of one-dimensional interacting fermion 
systems ("Luttinger liquid"), and one finds that the low-energy physical proper- 
ties are determined by only three parameters: the velocities of collective charge- and 
spin-density oscillations ("Up.o-) , and a coefficient Kp that determines the long-distance 
decay of correlation functions. These coefficients play a role similar to the Landau 
parameters of (three-dimensional) Fermi liquid theory. A number of physical proper- 
ties depending on these parameters are discussed below, but let us mention here that 
in particular the coefficient Kp is important in a much wider variety of phenomena: 
the temperature dependence of the NMR relaxation rate |0] or of X-ray scattering 
intensities the energy dependence of photoemission spectra 0, the effect of im- 
purities , or possible low-temperature ordered states in systems of coupled chains 
TT[| all depend on it. A brief discussion of bosonization will be given in chapter ^, 



but for more detailed and rigorous derivations and results, the reader is referred to 
more specialized articles [l^, |13|. 

As presented in sec. ^ the theory applies to translational invariant systems, 
and therefore one always has a conducting state. On the other hand, if one considers 
lattice systems, umklapp scattering breaks translational invariance. In sec. ^, I will 
show how this type of scattering can be incorporated into the bosonization formal- 
ism. One then finds insulating phases of the Mott-Hubbard type, with properties 
depending in an interesting way on the order of commensurability. There are two 
types of metal-insulator transitions: either at constant particle density, as a function 
of interaction strength, or at constant interaction, as a function of particle density. 
The "critical behavior" of the two types of transitions is quite different. 

A rather different approach (at least until recently) is based on the famous 



"Bethe ansatz" which in particular has made possible an exact solution of both 



continuum fermions interacting via (5-function potentials [^, |TB| and of the one- 
dimensional Hubbard model (and of many other interesting models). In section 
^ I will discuss the low-lying excitations of the one-dimensional Hubbard model as 
obtained from the exact solution. This will give a rather concrete illustration of the 
concept of "holons" and "spinous". Moreover, one observes interesting changes as 
the metal-insulator transition is approached. 

The exact Bethe ansatz eigenfunctions are so complicated that the direct cal- 
culation of correlation functions and many other physical properties of the one- 
dimensional Hubbard model is difficult even for very small systems [0 and impossible 
in the thermodynamic limit. In section O I present a method \TS\ that allows in 
particular a determination of the coefficient Kp for arbitrary correlation strength. 



2 



One then obtains a rather detailed and exact description of the low-energy (and low 
temperature) properties and also of the metal-insulator transition occurring when 
the average particle number per site, n, approaches unity. The method generalizes 
rather straightforwardly to other models, where however in general one has to rely 
on exact calculations for small systems and an extrapolation to the thermodynamic 
limit. 



2 Luttinger liquids 

2.1 The spinless case 

2.1.1 The model and its solution 

Let us consider the standard case of an interacting fermion model with the Hamilto- 
nian consisting of two parts: the kinetic energy and the interaction. Omitting spin 
for the moment, the kinetic energy term is of the form 

Hq ^^Skclck , (2.1) 

k 

where and c\ are the standard annihilation and creation operators for a fermion 
with momentum A;, and is the single-particle bandstructure. In a simple tight- 
binding model one would have — —2tcosk (the lattice constant is set to unity), 
but the precise form of Sk is unimportant here. The Fermi surface consists just of the 
two points i/cj?. 

For weak interactions between the particles, only states in the immediate vicin- 
ity of the Fermi points are important. For these states, one then can linearize the 
electronic dispersion relation around the Fermi points, and the kinetic energy term 
takes the form 

^0 = vp ^{(fc - kp)ala, + {-k - kF)blh} . (2.2) 

k 

Here the a (b) operators refer to states in the vicinity of +kF (— ^f), i-e. the a- 
particles move to the right, the 5-particles move to the left. The /c-summation is 
limited to an interval [— /cq, ^o] around kp (typically, ko ~ 7r/2, but the precise value 
isn't important here). The Fermi velocity is given by 



dsk 



(2.3) 

kp 



and the density of states is N{Ef) = 1/ {tivf)- In the Luttinger model, one generalizes 
this kinetic energy by letting the cutoff ko tend to infinity. There then are two 
branches of particles, "right movers" and "left movers", both with unconstrained 
momentum and energy. At least for weak interaction, this addition of extra states 
far from the Fermi energy is not expected to change the physics much. However, 
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this modification makes the model exactly solvable even in the presence of nontrivial 
interactions. Moreover, and most importantly, many of the features of this model 
carry over even to strongly interacting fermions on a lattice. 

We now introduce interactions between the fermions. As long as only forward 
scattering of the type {kp; —kp) — >• {kp; —kp) or {kp;kp) — >■ {kp;kp) is introduced, 
the model remains exactly solvable. Introducing the Fourier components of the par- 
ticle density operator for right and left movers by 

P+(?) = Il4+g«fc , P-(?)=II4+A ' (2-4) 

k k 

the interaction Hamiltonian describing these processes takes the form 

H^nt = ^T.{^92{q)pAq)p~i-Q) + 9dq)[p+iq)p+i-q) + p-i-q)p-iq)]} -(2.5) 

Here, g2iq) and g4,{q) are the Fourier transforms of a real space interaction potential, 
and in a realistic case one would of course have g2{q) = g^iq), but it is useful to allow 
for differences between g2 and (74. For Coulomb interactions one expects (72, ^74 > 0. 
In principle, the long-range part of the Coulomb repulsion leads to a singular q- 
dependence. Such singularities in the gi can be handled rather straightforwardly and 



can lead to interesting physical effects [19|, but here I shall limit myself mainly to 



nonsingular 5'2;5'4- Electron-phonon interactions can lead to effectively attractive 
interactions between electrons, and therefore in the following I will not make any 
restrictive assumptions about the sign of the constants. One should however notice 
that a proper treatment of the phonon dynamics and of the resulting retardation 
effects requires more care pO| . 



The model defined by eqs. ( |2.2| ) and ( |2.5| ) can be solved exactly. The solution 
is based on the following facts: 

1. the density fluctuation operators p± obey Bose type commutation relations: 

[p4-q),p4q')] = [p-(?), p- (-?')] = v|§ , [p+(?), p- (?')] = o .(2.6) 

Moreover, for g > both p+(— g) and P-{q) annihilate the noninteracting 
groundstate. These properties are closely related to the existence of an in- 
finity of states at negative energies and would not have been true in a lattice 
system like the one described by (|2.1|). 



2. The kinetic part of the Hamiltonian can be re-written as a term bilinear in 
boson operators, i.e. quartic in fermion operators: 

^0 = ^T.[pMp^i-q) + p-i-q)p-iq)] ■ (2.7) 
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This equivalence may be made plausible noting that p+ (g) creates particle-hole 
pairs that all have total momentum q. Their energy is ek+q — £k, which, because 
of the linearity of the spectrum equals vpq, independently of k. Thus, states 
created by p+{q) are linear combinations of individual electron-hole excitations 
all with the same energy, and therefore are also eigenstates of ( p.2|) . 

3. The states created by repeated application of p± on the ground state form a 
complete set of basis states [0, 

Putting together ( |2.7|) and ( |2.5| ), the complete interacting Hamiltonian then becomes 
a bilinear form in boson operators, that is easily diagonalized by a Bogolyubov trans- 
formation. A first consequence is the expression for the excitation spectrum 

^(g) = MivF + 9,{q) /{27r)f - {g,{q) / i2n)YY/' . (2.8) 

The diagonal boson operators are linear combinations of the original p operators, and 
consequently, these elementary excitations are collective density oscillations, their 
energy being determined both by the kinetic energy term and the interactions. 

We note here that in order for the Bogolyubov transformation to be a well- 
defined unitary transformation, g2{q) has to decrease at large q at least as 
on the other hand, the large-g behavior of g2 is unimportant for the low-energy 
properties of the model. We therefore in the following will almost always use a 
g-independent g2 and g^. An approximate and frequently used way to cure the 
divergences arising due to this procedure is to keep the parameter a in subsequent 
formulae as a finite short-distance cutoff, of the order of a lattice spacing. One 
can then also include the "backward scattering" {kp; —kp) — > {—kp', kp), because for 
spinless electron this is just the exchange analogue of forward scattering and does 
not constitute a new type of interaction. 

Up to this point, this construction does not allow for a direct calculation of 
correlation functions like the one-particle Green's function or more generally any 
function involving individual creation or destruction operators. This type of correla- 
tion function becomes tractable by representing single particle operators in terms of 
the boson operators; To this end, we introduce the field operators 

m = -TE-e'"'''^'"^'^[P+(p)+P-(p)]-^^ ' (2-9) 
^ P^oP ^ 

n(x) = jJ2e~'^^P^/'-^^^[pM-P-ip)] + J/L . (2.10) 

Here and J are the number of particles added to the ground state and the difference 
between the number of right and left-moving particles, respectively, and a is a cutoff 
parameter which (at least in principle, see the discussion above) has to be set to zero 
in the end of any calculation. and 11 obey canonical boson commutation, relations: 

[ct>{x),Il{y)]=t6{x-y) , (2.11) 
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and is related to the local particle density via = —irp^x). The expression 

for the single fermion operators then is 

1 



'0±(a:) 



lim 



: exp [±ifci?x =F i(j){x) — i6{x)] 



(2.12) 



-0 y/27ia 

where 6{x) = it Il{x')dx', and the upper and lower sign refer to right- and left- 
moving electrons, respectively. A detailed derivation of this important relation as an 



operator identity is given in the literature |]T3|, However, a simple plausibility ar- 
gument can be given: creating a particle at site x means introducing a kink of height 
71 in 0, i.e. at points on the left of x (p has to be shifted by vr. Displacement oper- 
ators are exponentials of momentum operators, and therefore a first guess would be 
ip'^i^x) ~ exp(i7r J^^Il{x')dx'). However, this operator commutes with itself, instead 
of satisfying canonical anticommutation relations. Anticommutation is achieved by 
multiplying with an operator, acting at site x, that changes sign each time a particle 
passes through x. Such an operator is exp(±20(x)). The product of these two factors 
then produces ( p.29| ). 

The full Hamiltonian can also be simply expressed in terms of (p and H. Ne- 
glecting the momentum dependence of the gi, one easily finds 



H = Hq + Hj 



int 



dx 



-mxY + 



(2.13) 



2 ' ' 27cK 

This is obviously just the Hamiltonian of an elastic string, with the eigenmodes corre- 
sponding to the collective density fluctuations of the fermion liquid. It is important to 
notice that these collective modes are the only (low-energy) excited states, and that 
in particular there are no well-defined single particle excitations. The parameters in 
( p. 13 ) are given by 



u 



l{vF + g,/{27^)r-gl/{27rn/' , K 



2tivf + g4- g2 
2ttvf + gi + g2 



1/2 



(2.14) 



The energies of the eigenstates are uj{q) = u\q\, in agreement with eq. (|2.8|). 



2.1.2 Physical properties 

The simple form of the Hamiltonian ( ^.13| ) makes the calculation of physical proper- 
ties rather straightforward. First note that acoustic phonons in one dimension have 
a linear specific heat. Consequently, the low-temperature specific heat of interacting 
fermions is C (T) = 7T, with 

7/70 = vf/u . (2.15) 

Here 70 is the specific heat coefficient of noninteracting electrons of Fermi velocity 
vf- Also, the coefficient u/K fixes the energy needed to change the particle density, 
and consequently the renormalization of the compressibility k, is given by 

k/kq = vfK/u , (2.16) 
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where kq is the compressibihty of the noninteracting case. 

The quantity 11 is proportional to the current density. Obviously, the Hamil- 
tonian commutes with the total current, and therefore the frequency dependent con- 
ductivity is a delta function at c<j = 0. Using the Kubo formula, one straightforwardly 
finds 

a{uj) = Ku5{uj) , (2.17) 

i.e. the product Ku determines the weight of the dc peak in the conductivity. 

The above properties are those of an ordinary Fermi liquid, the coefficients m, 
and K determining renormalizations with respect to noninteracting quantities. We 
will now consider quantities which show that a one-dimensional interacting fermion 
system is not a Fermi liquid. Consider the single-particle Green's function which can 
be calculated using the representation (|2.12|) of fermion operators: 

G«(x,t) = -?^(t)([^+(a;,t), 4(0,0)]+) 

.2 1<5/2^ 



e^'^^^Re 



TT \ Ut — X 



a 



(a + lutY + x"^ 



(2.18) 



One then finds for the momentum distribution function in the vicinity of fci?: 

Uk ^ — consi.sig'n.{k — kF)\k — kF\^ , (2-19) 

and for the single-particle density of states: N{uj) ~ |a;|'^, with 5 = {K + 1/K — 2)/4, 
and /9 is a model-dependent constant. Note that for any K 1, i.e. for any 
nonvanishing interaction, the momentum distribution function and the density of 
states have power-law singularities at the Fermi level, with a vanishing single particle 
density of states at Ep. This behavior is obviously quite different from a standard 
Fermi liquid which would have a finite density of states and a step-like singularity 
in n^. The absence of a step at kp in the momentum distribution function implies 
the absence of a quasiparticle pole in the one-particle Green's function. Instead one 
finds in the spectral function A{k, uj) a continuum above a threshold u{k — kp), with 
a singularity close to the threshold: A{k,Lj) oc (cj — u{k — kp)Y~^, and there is also 



some (non-divergent) spectral weight at negative energies p^ , |23| , |24|. 

The coefficient K also determines the long-distance decay of all other correlation 
functions of the system: in the present context the most interesting correlations are 
those involving charge density or pairing fiuctuations. Using the representation ( p. 121 ) 
the corresponding operators can be written as 

OcDwix) = V'l(a:)V'+(a;) = lim^^e-2^'^(^) , (2.20) 
Osc{x) = ^/^_(x)7/'+(x) = lim— e-2*^(^') . (2.21) 
Similar relations can also be found for other operators. The CDW and pairing corre- 



lation functions now are easily calculated as a function of x, t, and temperature [12 . 
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I here only give the leading asymptotic behavior at zero temperature: 

{Olj,w{x,t)OcDwm = Ciix^-uHY'' , 

{Olc{x,t)Oscm = C^ix'-uH'r'/^ , (2.22) 

with interaction-dependent constants Ci,2- The corresponding susceptibilities (i.e. 
the Fourier transforms of the above correlation functions) behave at low temperatures 
as 

XcDwiT) ^ T2^-2 , xsciT) ^ T2/^'-2 , (2.23) 

i.e. for K < 1 charge density fluctuations at 2^^;' are enhanced and diverge at low 
temperatures, whereas for K > 1 pairing fluctuations dominate. The remarkable fact 
in all the above results is that there is only one coefficient, K, which determines all 
the asymptotic power laws. 

To summarize the spinless case, we have found striking differences with the usual 
Fermi liquid picture: there are no quasiparticle poles in the single-particle Green's 
function, but rather power law singularities, with interaction-dependent exponents. 
Similar power laws appear in all types of correlation functions. However, all exponents 
are related to one parameter K, which contains all the interaction dependence. This 
type of behavior has been called Luttinger liquid by Haldane [|13] . 



2.2 Spin— 1/2 fermions 

In the case of spin-1/2 fermions, all the fermion operators acquire an additional spin 
index s. Following the same logic as above, the kinetic energy then takes the form 

Ho = Vf ^{{k - kF)a\ sa'k,s + {-k - kF)h\ j)k,s] 



k,s 

2ttvf 



^ [p+Mp+A-(D + P-A-(l)P-A(D] , (2.24) 

^ -?>0,s 

where density operators for spin projections s =], | have been introduced: 

P+,s(g) = I]aI+g,.aM . P-A^l) = H^l+q,sKs ■ (2.25) 

k k 

There are now two types of interaction. First, the "backward scattering" [kF-, s; — fc^, t) 
—>■ {—kF, s; kF,t) which for s ^ t cannot be re-written as an effective forward scat- 
tering (contrary to the spinless case). The corresponding Hamiltonian is 

Hint,l = ^ X! 9lO'k,sK,t^P+'^kp+q,tbk-2kp-q,s ■ (2.26) 
kpqst 

And, of course, there is also the forward scattering, of a form similar to the spinless 
case 

Hrnt,2 = -^E{2^?2(g)p+,s(g)p-,t(-g) (2.27) 

qst 

+g4{q)[p+A(l)p+A~Q)+P-A-'l)p~A'l)]} ■ (2-28) 
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To go to the bosonic description, one introduces and 11 fields for tlie two spin 
projections separately, and then transforms to charge and spin bosons via (pp^a = 
(01 ± (j)i)/y/2, IIp^o- = (Hi ± The operators (pi, and Hi, obey Bose-like 

commutation relations: 

[(puix), Upiy)] = i6uf,6{x - y) , 
and single fermion operators can be written in a form analogous to (|2.12|) : 



lim 



V2 



: exp 



±ikFX-ii±i(l)p + s(l)^)-i9p + s9^))/V2\ , (2.29) 



where 9u{x) = vr Il,^{x')dx' . 

The full Hamiltonian H = Hq + Hint i + Hint 2 then takes the form 



H = Hn + H„ 



2gi 



dx cos 



(27ra)2 

Here a is a short-distance cutoff, and for u 



p,a 



dx 



with 



27rvF + 2gf4^^ + gf^ 
27CVF + 25(4 . - g^, 



1/2 



(2.30) 



(2.31) 



(2.32) 



and gp = gi — 2g2, g^ = gi, gi,p = g^, g4,a = 0. For a noninteracting system one 
thus has Ui, = vf (charge and spin velocities equal!) and K^, = 1. For gi = 0, ( p.30D 
describes independent long-wavelength oscillations of the charge and spin density, 
with linear dispersion relation ujy{k) = u^\k\, and the system is conducting. As in 
the spinless case, there are no single-particle or single particle-hole pair excited states. 
This model (no backscattering) , usually called the Tomonaga-Luttinger model, is the 
one to which the bosonization method was originally applied [E^, [ 



For gi ^ the cosine term has to be treated perturbatively. A straightforward 
calculation gives the renormalization group equation 



4i(o = ^^7i(o^ ,42(0 

dl txvf dl 



-gM 



(2.33) 



21TVF 

where the renormalized cutoff Ec is related to the bare cutoff Eq via Ec = Eoexp{l). 
Thus, for repulsive interactions {gi > 0), gi is renormalized to zero in the long- 
wavelength limit, and at the fixed point one has K* = 1. The three remaining 
parameters in ( ^.3U ) then completely determine the long-distance and low-energy 
properties of the system. 



It should be emphasized that ( |2.30D has been derived here for fermions with 
linear energy-momentum relation. For more general (e.g. lattice) models, there 
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are additional operators arising from band curvature and the absence of high-energy 
single-particle states. One can however show that all these effects are, at least for not 
very strong interaction, irrelevant in the renormalization group sense, i.e. they do not 
affect the low-energy physics. Thus, l \2. 3Q ) is still the correct effective Hamiltonian 
for low-energy excitations. 

The Hamiltonian ( p. 301) also provides an explanation for the physics of the case 
of negative gi, where the renormalization group scales to strong coupling (eg. ( p. 33]) ). 
In fact, if \gi\ is large in ( p.30|) , it is quite clear that the elementary excitations of 
(per will be small oscillations around one of the minima of the cos term, or possibly 
soliton-like objects where (pa goes from one of the minima to the other. Both types of 
excitations have a gap, i.e. for gi < one has a gap in the spin excitation spectrum, 
whereas the charge excitations remain massless. This can actually investigated in 



more detail in an exactly solvable case |28 



2.2.1 Spin— charge separation 

One of the more spectacular consequences of the Hamiltonian ( p.30| ) is the complete 
separation of the dynamics of the spin and charge degrees of freedom. For example, in 
general one has Uo- 7^ Up, i.e. the charge and spin oscillations propagate with different 
velocities. Only in a noninteracting system or if some accidental degeneracy occurs 
does one have = Up = vp- To make the meaning of this fact more transparent, let 
us create an extra particle in the ground state, at t = and spatial coordinate xq. 
The charge and spin densities then are easily found, using p{x) = —{\/2/Ti)d(pp/dx 
(note that p{x) is the deviation of the density from its average value) and az{x) = 
-{V2/7r)d^Jdx : 

{0\ip+{xo)pix)ip\.{xo)\0) = 5{x-xo) , 

(0|V'+(xo)a,(x)4(xo)|0) = Six-xo) . (2.34) 

Now, consider the time development of the charge and spin distributions. The time- 
dependence of the charge and spin density operators is easily obtained from ( p.30| ) 
(using the fixed point value gi =0), and one obtains 

{0\^+{xo)p{x,t)^l{xo)\0) = 5{x - Xo-Upt) , 

(0|^/'+(xo)(T2(x,t)V'|(xo)|0) = 5{x-Xo-Uat) . (2.35) 

Because in general 7^ Up, after some time charge and spin will be localized at 
completely different points in space, i.e. charge and spin have separated completely. 
A interpretation of this surprising phenomenon in terms of the Hubbard model will 
be given in sec.(^. 

Here a linear energy-momentum relation has been assumed for the electrons, 
and consequently the shape of the charge and spin distributions is time-independent. 
If the energy-momentum relation has some curvature (as is necessarily the case in 
lattice systems) the distributions will widen with time. However this widening is 
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proportional to -y/t, and therefore much smaller than the distance between charge 
and spin. Thus, the qualitative picture of spin-charge separation is unchanged. 



2.2.2 Physical properties 

The simple form of the Hamiltonian ( ^.3(J| ) at the fixed point gl = makes the 
calculation of physical properties rather straightforward. The specific heat now is 
determined both by the charge and spin modes, and consequently the specific heat 
coefficient 7 is given by 



7/70 = T^i^p/Up + Vp/U^) 



(2.36) 



Here 70 is the specific heat coefficient of noninteracting electrons of Fermi velocity 

The spin susceptibility and the compressibility are equally easy to obtain. Note 
that in ( |2.3U| ) the coefficient UajK^ determines the energy necessary to create a 
nonzero spin polarization, and, as in the spinless case, Up/Kp fixes the energy needed 
to change the particle density. Given the fixed point value K* = 1, one finds 



x/Xo = vf/u^ , k/kq = vfKp/u 



p ' 



(2.37) 



where xo and are the susceptibility and compressibility of the noninteracting case. 
From eqs.( p.36| ) and (|2.37| ) the Wilson ratio is 



p _ A ;u 
riy/ — 



X To 
7X0 



Up + u„ 



(2.38) 



The quantity np(x) is proportional to the current density. As before, the Hamil- 
tonian commutes with the total current, one thus has 



g{uS) = 2KpUpb{uS) + cr. 



reg 



(2.39) 



i.e. the product KpUp determines the weight of the dc peak in the conductivity. The 
regular part of the conductivity in general varies as ur' at low frequencies p9| . 

The above properties, linear specific heat, finite spin susceptibility, and dc con- 
ductivity are those of an ordinary Fermi liquid, the coefficients Up,Ua, and Kp de- 
termining renormalizations with respect to noninteracting quantities. However, the 
present system is not a Fermi liquid. This is in fact already obvious from the preceding 
discussion on charge-spin separation, and can be made more precise considering the 
single-particle Green's function. Using the representation ( |2.29|) of fermion operators 
one finds (at the fixed point gi = 0) 



G^(x,t) = -^^(t)([V^+,s(a:,t),4,.(0,0)]+) 



9{t) 



Re 



TT 



'{Upt — x){Uat — x) 



a 



[a 



+ iupty 



5/2" 



(2.40) 
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Note that this expression factorizes into a spin and a charge contribution which 
propagate with different velocities. Fourier transforming ( p. 40 ) gives the momentum 
distribution function in the vicinity of kp: 

Uk ^ — const. sign{k — kF)\k — kpl^ , (2-41) 

and for the single-particle density of states (i.e. the momentum-integrated spectral 
density) one finds: 

N{uj) ^ \uj\^ . (2.42) 

In both cases 6 = {Kp + l/Kp — 2)/4. Note that for any Kp ^ 1, i.e. for any 
nonvanishing interaction, the momentum distribution function and the density of 
states have power-law singularities at the Fermi level, with a vanishing single particle 
density of states at Ep. This behavior is obviously quite different from a standard 
Fermi liquid which would have a finite density of states and a step-like singularity 
in nfc. The absence of a step at fci? in the momentum distribution function implies 
the absence of a quasiparticle pole in the one-particle Green's function. In fact, a 
direct calculation of the spectral function A{k, u) from ( |2.4CI| ) p3| , [2^ shows that 
the usual quasiparticle pole is replaced by a continuum, with a lower threshold at 
min(Mj,)(fc — kp) and branch cut singularities aX oj = UpP and u = u„p: 

A{k,u})^ {uj-u^{k-kp)Y'^/^ ,\uj-Up{k-kp)\'^^-^'>/^ {up>u^) , (2.43) 

A{k,u})^{uj-Up{k-kpyf-^^/Wu}-u^{k-kp)f-^/^ {up<u^) . (2.44) 

The coefficient Kp also determines the long-distance decay of all other corre- 
lation functions of the system: Using the representation ( |2.29| ) the charge and spin 
density operators ai 2kp are 

2ikpx 

OcDwix) = ^^l,,(x)V^+,,(x) = lim e-*^<^''(^)cos[y20^(x)] , (2.45) 

^ a^o Tea 

OsDwM = E^^s(^)^+,-(^) = l™^^e-^"^*''(^^cos[V2e.(x)] . (2.46) 

Similar relations are also found for other operators. It is important to note here 
that all these operators decompose into a product of one factor depending on the 
charge variable only by another factor depending only on the spin field. Using the 
Hamiltonian ( p.30| ) at the fixed point gl = one finds for example for the charge and 
spin correlation functions^ 

{n{x)n{0)) = Kp/{7rxy + Aicos{2kpx)x-^-^''\n~^^^{x) 

+ A2Cos{Akpx)x-^^^ + ... , (2.47) 
(5(x)-5(0)) = l/{7Txf + Bicos{2kpx)x-^^^'ln^/\x) + ... , (2.48) 



*The time- and temperature dependence is also easily obtained, see E 
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with model dependent constants Ai,Bi. The elhpses in ( |2.47] ) and ( |2.48| ) indicate 
higher harmonics of cos(2A;ira;) which are present but decay faster than the terms 
shown here. Similarly, correlation functions for singlet (SS) and triplet (TS) super- 
conducting pairing are 

{0!rsS^)OTsAO)) = Dx-'-'/^^ln'/\x) + ... . (2.49) 
The logarithmic corrections in these functions ||30| have been studied in detail recently 



|3l| , |32| , p^ , |3^ . The corresponding susceptibilities (i.e. the Fourier transforms of the 



above correlation functions) behave at low temperatures as 

XcDiy(T)^T^''-i|ln(T)r3/2 , x5w(T) ^ T^'^-^l ln(T)| , (2.50) 
X55(T)^TV^''-^|ln(T)r3/2 ,;^^^(T)^Ti/^p-i|ln(T)|V2 , (2.5I) 

i.e. for Kp < 1 (spin or charge) density fluctuations at 2kp are enhanced and diverge 
at low temperatures, whereas for Kp > 1 pairing fluctuations dominate. These cor- 
relation functions with their power law variations actually determine experimentally 
accessible quantities: the 2k p and 4k f charge correlations lead to X-ray scattering 
intensities l2kp ~ T^f", I^kp ~ y^Xp-i^ similarly the NMR relaxation rate due 
to 2kF spin fluctuations varies as 1/Ti ^ T^p. The remarkable fact in all the above 
results is that there is only one coefficient, Kp, which determines all the asymptotic 
power laws. 

Correlation functions in the spin-1/2 case ("spin-1/2 Luttinger liquid") share 
one important property with spinless fermions: they have power-law behavior, with 
interaction-dependent powers determined by one coefficient, Kp. However, the phe- 
nomenon of spin-charge separation adds some additional features in this case and 
has spectacular consequences both for thermodynamical and spectral properties. 



3 Umklapp scattering and metal— insulator transi- 
tions 

3.1 Half-filling 

In the model discussed in the preceding section, total momentum was always a con- 
served quantity, and consequently these models are metallic with an infinite dc con- 
ductivity. However, in a half-filled band one has kp = vr/2. Then umklapp scatter- 
ing, transferring two particles from —kp to kf, involves momentum transfer Akp = 27r 
which is a reciprocal lattice vector. These processes are thus allowed and lead to an 
extra term 

Hint,3 = ^ X! 93{(^k,s^p-sbp-'ikF+q-sbk-2kF-q,s + h.C.) . (3.1) 
kpqs 
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in the Hamiltonian. Note that because of the Pauh principle, only scattering of 
two particles of opposite spin is allowed if, as usually is assumed, is momentum- 
independent. In the boson representation, this term leads to an additional interaction 
for the charge degrees of freedom: 

Hint,?. = TTT^ / dxcos{V8(j)p) . (3.2) 

Similarly to the gi term, this term can be handled via a renormalization group cal- 
culation. The equations are 

^93{l) = -^(giil) - ^g^g^il) , ^.{gi{l) - 2^2(0) = -^Q^ilf ■ (3-3) 

These in fact are the well-known Kosterlitz-Thouless equations. There are two 
regimes: (i) for Ig^l < gi ~ 2g2 gs scales to zero, and gi — 2g2 to the fixed point 
value {gi — 2g2)* = [{gi — 2g2Y — gf^^^"^- Consequently, the charge excitation spec- 
trum remains massless, but there are corrections of order g^ to Kp-. in eq ( p.32| ) for 
Kp one replaces gi — 2g2 —>■ {gi — 2g2)*, but Up is unrenormalized at this order. Of 
course, in this case the system remains a metal, (ii) for {g^l > gi — 2g2 gs scales 
towards strong coupling. From eq. (|3^) one then expects 0p to be essential fixed to 
one of the minima of the cosine potential. This gives rise to a gap in the charge 
excitation spectrum |^ ||. The ground state then is insulating. 



In a particular but illuminating case, the umklapp problem can be solved ex- 
actly [^: after the unitary transformation 0p — 0p/-\/2, Up V^IIp, the total 
Hamiltonian for the charge degrees of freedom becomes 

Hp = jdx nUpKpUl + £^(40p)2 + I dxcos{2^p) . (3.4) 

Using the transformations discussed in sec. 0, this problem can be transformed into 
a spinless fermion model. In particular, for Kp = 1/2, the corresponding spinless 
model has K = 1 (cf. eq. (p.l3|) ), i.e. g2 = 0. Using eq.( p.20| ) Hp then transforms 
into 

Hp = VFy2{{k-kF)alak + {-k-kF)blbk} + ^^y2{albk^2kp + blak+2kp) -(3.5) 
k 27ra ^ 

This form is easily diagonalized, and one finds an excitation spectrum for the spinless 
fermions of the form (noting that for half-filling kF = v/2) 

Efc = ±[4(A;±7r/2)2 + AY/' , (3.6) 
with a gap 

A = W(2vra) • (3-7) 

In the ground state all negative energy states are filled, all positive energy states 
are empty, and because of the gap the system then is an insulator. Note that the 
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fermion operators in ( p.5| ) in fact create solitons in the 0p field, without affecting (p^. 
At Kp 7^ 1/2, the physical picture is not changed qualitatively, in particular there is 
still a gap in the charge excitation spectrum, however, the functional dependence is 
changed: in general, one has A oc gi^, with u = 1/(2 — 2Kp). 

The spin field 0^ is completely unaffected by the metal-insulator transition. 
Consequently, even the insulating state has a Pauli-like susceptibility, and long- 



range spin correlations of the type (|2.48|) , with Kp = 0. This behavior is that of 



the one-dimensional antiferromagnetic Heisenberg model, which can be considered 
of electrons localized on individual atoms. In the present case, the gap is 
typically much smaller than the bandwidth, and the localization length is therefore 
rather big. Nevertheless, the qualitative behavior of the correlations is the same. 

We thus here have an example of a metal-insulator transition, occuring at 
constant particle density (half-filling) as a function of the interaction parameters. 
Approaching the transition from the metallic side, one has Kp ^ 1, Up ^ const., 
and thus in particular the Drude weight of the conductivity remains finite and jumps 
discontinuously to zero as the transition line is crossed. This jump is nothing but 
the usual jump of the stiffness parameter (or "superfluid density") of the Kosterlitz- 
Thouless transition. On the insulating side, the gap opens exponentially as a function 
of interaction strength, again the typical Kosterlitz-Thouless behavior. 

One clearly would like to understand what happens if one is close to but not 
exactly half-filled. One thus adds a chemical potential term to the original model to 
obtain 

3 

H = Ho + Y^ Hint,i - lJ^J2("'i,s(^k,s + bi^sKs) ■ (3.8) 

i=l k,s 

In the bosonic description, the chemical potential term gives rise to a term pro- 
portional to d(j)p/dx, which in turn, in the spinless fermion language of eq.( p.5|) is 
equivalent to a chemical potential. For Kp = 1/2, the solution is still straightforward: 
as long as < A, the ground state is unchanged, but for /i > A, positive energy 
states start to be occupied (and similarly, for yu < —A, negative energy states are 
emptied). The number of extra carriers (i.e. the deviation of the carrier density from 
half-filling) varies as 6n oc — A. The system is now metallic, and from eq.( |3.6| ) 
one obtains the effective velocity of the charge modes as 

^ vl\kF-7c/2\ 
' [4(fc^-7r/2)2 + A2]V2 ^'"•^^ 

In particular, Up vanishes linearly as half-filling is approaches. Thus, the Drude 
weight of the conductivity vanishes linearly as n — 1. 

If Kp 7^ 1/2 in eq.( p.4| ) there are additional interaction terms between the 
spinless fermions of eq. ( ^.51 ). However, in the vicinity of half-filling, these interactions 



can be eliminated from the problem |]36| , p7[ , and one has then for the effective 
parameter governing the low-energy physics K* = 1/2 in all cases, even if the original 
Kp was quite different from 1/2. 
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It should be emphasized here that K* = 1/2 is valid for n — 1, but not for n = 1. 
In the latter case, there is a gap in the charge excitation spectrum due to the umklapp 
term {\i.'2\), and the correlations of 0p become long ranged, i.e. K* = 0. Close to 
half-filling, the asymptotic behavior of the charge part of correlation functions like 
( p.48|) is essentially determined by the motion of the added carriers. Writing the 
density of carriers as p = \1 — n\ one then expects a crossover of the form 

■ 5(0)) ^ cos(2A;^x)[l + (px)^]-^''/^^-! \n^/\x) (3.10) 

for the 2kF part of the spin correlation function, and similarly for other correlation 
functions. Clearly, only for x ^ 1/p are the asymptotic power laws valid, whereas 
at intermediate distances 1 <^ x ^ 1/p one has effectively Kp = 0. The form ( p.lO] ) 
provides a smooth crossover as n 1. 

An interesting question is the sign of the charge carriers, especially close to the 
metal-insulator transition. The standard way to determine this, the sign of the Hall 
constant, is useless in a one-dimensional system. As an alternative, the thermopower 
can be used which is negative (positive) for electron (hole) conduction. In general, 
calculation of the thermopower is a nontrivial task, as the curvature of the bands plays 
an important role, and the approximate form of the Hamiltonian ( |2.30| ) is therefore 
insufficient. Moreover, both charge and spin entropies can play a role. However, close 
to the metal-insulator transition Up -C Mo-, and therefore the entropy of the charge 
degrees of freedom is much bigger than the spin entropy. As discussed above, the 
charge part of the Hamiltonian can be transformed into a model of massive fermions, 
with energy-momentum relation given by (|3.6| ). At half-filling all negative energy 
states are filled, all positive energy states are empty. Doping with a concentration 
n* of holes, some of the negative energy states become empty and only states with 
\k\ > kp (X n* are filled. Because of the vanishing interaction, a standard formula for 



the thermopower can be used |38| and gives 



~ 6|e| t;2(A;^)2(t;2(A;^)2 + A2)l/2 ' (3.11) 

i.e. approaching the metal-insulator transition from n < 1, the thermopower is hole- 
like, whereas obviously far from the transition (ra <^ 1) it is electron-like. The exactly 
opposite behavior can be found for n > 1. 

At zero temperature and away from half-filling, one has an infinite dc conduc- 
tivity in this model. However, at finite temperature there is some probability to 
excite a carrier into a momentum state that makes umklapp scattering possible. One 
then expects a conductivity increasing exponentially at low temperature p^ j. 

For spinless fermions, a term like ( p.l|) cannot act, because of the Pauli principle. 
However, a term like 

dx {ipUx)[d^ipl{x)]ip.{x)[d^ip.{x)] + h.c.} ^ J dx cos(40(x)) (3.12) 

does produce umklapp scattering at half-filling and can give rise to a gap in the 
charge excitation spectrum if i^' < 1/4 ||39|, The metal-insulator transition 
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has properties analogous to that in the spin-1/2 case and has been studied in detail by 
Shankar |42|. A lattice model having this type of transition is the XXZ spin chain, 



which via a Jordan-Wigner transformation can be seen as a model of interacting 
spinless fermions. 



3.2 Other commensurabilities 

The umklapp operator ( |3.2| ) is effective at and near half-filling. One might wonder 
whether higher-order "commensurability" , i.e. in a third- or quarter-filled band, can 
also lead to insulating states, and at least intuitively it seems clear that this should 
be possible. On the other hand, for such cases, simple two-particle processes cannot 
be simultaneously momentum conserving (even modulo a reciprocal lattice vector) 
and only involve states near the Fermi energy. However, many-particle processes 
induced by bare two-particle interactions do exist: consider the case of a third-filled 
band {kp = 7r/3). In a first step then one has a (non-umklapp) process of the type 
(— 7r/3; — 7r/3) (— 7r;7r/3), with the particle at — tt in a high-energy state at the 
band edge. In a second step one makes the transition (— tt = tt; — 7r/3) — * (7r/3; 7r/3). 
Umklapp here intervenes due to the identification of states at — tt and tt, and the 
final result is the transfer of three particles from —kp to kp. For weak two-particle 
interactions with an interaction matrix element g the effective matrix element Qu for 
this type of process is of order g'^ /t (with bandwidth t), but for stronger interactions 
it is hard to estimate. The generalization to bandfilling 1/m (i.e. 2/m particles 
per site) is straightforward: one transfers m particles from —kp to kp, with matrix 
element ~ g"^^^ /t"^. 

To take the existence of this type of processes into account in the original 
continuum description, one adds a term 

H^ = gu j dx{[i,l{x)r[^^{x)r + h.c.} , (3.13) 

where the product of many fermion operators at one site is to be understood as 
point-split, and the spin summation is implied. In the bosonic language, this term 
translates into 

= gujdx{OcDw{xT + h.c.} (3.14) 

^ gu J dx cos"\V2(p^{x)) cos{V2m(j)p{x)) . (3.15) 

There are now two physically rather different cases, according to whether n is even or 
odd. Consider first the even case. Then the cos"^{\/2(f)f^{x)) term can be expanded, 
and the lowest-order (most relevant in the renormalization group sense) comes from 
the constant in this expansion, e.g. the effective Hu is 

Hu gu dx cos{V2m(f)p{x)) . (3.16) 
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After a simple unitary transformation rescaling 0p, we have now a problem formally 
identical to the half-filled case. At filling exactly equal to 1/m we thus have an 
insulating or a metallic phase, with the metallic state stable for Kp > and 
K* = 4/m'^,u* = const, at the metal-insulator transition. For m > 2 the condition 
for the existence of an insulating state, Kp < usually requires strongly repulsive 

interactions, and even then can not always be satisfied: e.g. for the Hubbard model, 
Kp > 1/2 even for infinitely strong repulsion. In the insulating state, 0p has only 
small oscillations around a fixed average value, and thus the effective Kp is zero. We 
then have long-range charge density wave order at wavevector Akp, but still the spin- 
spin correlations typical of the one-dimensional antiferromagnetic Heisenberg model 
(cf. eq. (|2.48|) ). We note that the solitons in this case have a jump in 0p of A/27r/m, 
and thus carry charge 2/m. As in the half-filled case, the magnetic properties are 
largely independent on whether one is in the insulating or metallic state. 

Varying the particle density at constant Kp < one recovers behavior 

similar to the nearly-half-filled case above: as n ^ 2/m one has Up ^ and Kp — >■ 
2/m^. Of course, for m = 2 we recover precisely the results of the preceding section. 

Things are a bit more complicated if m is odd. In this case the most relevant 
term in Hu takes the form 



9u 



dx cos(-\/20o-(a;)) cos(-\/2m0p(x)) 



(3.17) 



e.g. spin and charge degrees of freedom are coupled. The scaling dimension of this 
operator is (1 + m'^Kp)/2, and consequently this operator is relevant and produces 
a gap if Kp < The insulating state is characterized by a large effective gu, 

and consequently, from ( p.l7| ) one expects both (pp and 0o- to oscillate around stable 
equilibrium positions, i.e. there is a gap in both the charge and the spin excitations, 
and the ground state is non-magnetic and one has 2kF CDW ordering (in contrast 
to the paramagnetic ground state and Akp CDW found for m even). There are now 
separate soliton-like excitations for charge and spin, carrying charge 2/m and spin 
1/2, respectively. 

In a particular case, an exact solution of the odd-m problem can be given: 
assume gi = (which anyway is the fixed-point value for repulsive interactions). 
Then, after the unitary transformation (pp (pp/m, the total Hamiltonian, e.g. the 
sum of ( |2.30| ) and ( |3.17| ) becomes 



H 



dx 
+ I dx 



TTUpKpm"^ 



2TrKpm'^ 



+gu / dx cos(v^0ct(x)) cos(V20p(x)) , 



(3.18) 



The 5f„-term now has exactly the form of the 2k p CDW operator ( |2.45| ) 
bilinear in fermion operators. Moreover, for Kp = l/m^ and Up = u^r 



i.e. it is 

= u, the 
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first two terms in ( p.l8| ) represent free spin-1/2 fermions, e.g. in tliis case tlie full 



Hamiltonian is in fact the spin-1/2 version of 

k,s 

+9uJ2((^lsh-2kp,s + bl,ak+2kp,s) ■ (3.19) 



k,s 



This form is of course easily diagonalized, and, as expected one has a gap both in the 
charge and in the spin excitations. In this particular case these gaps are equal, but 
generally, e.g. for Kp ^ this will not bee the case. 

If one is away from filling 1/m, an additional chemical potential term appears 
in ( p.l9| ), and either states above the gap get filled or states below the gap become 



empty. This is very similar to the nearly-half-filled case. Because we have effectively 
noninteracting electrons in ( |3.19| ), one still has Kp = 1/m^, and Up oc \n — 2/m\. 



However, and contrary to the case of m even, now the spin velocity shows critical 
behavior: oc \n — 2/m\, leading in particular to a diverging spin susceptibility. 

This behavior is however not generic: if Kp deviates only slightly from 1/m^, or 
one of the other solvability conditions is not satisfied, there are interaction terms in 



addition to the free fermion Hamiltonian (|3.19| ). Close to filling 1/m there are very 



few extra fermions, compared to n = 2/m. However, and contrary to the spinless 
case relevant for m even, interactions do have dramatic effects on spin-1/2 fermions 
even in the very dilute limit. In fact, as long as the interaction is not long-ranged, 
for sufficiently high dilution (more precisely, if the interparticle distance is bigger 
than the scattering length), the precise form of the interaction matrix element is 
expected to be unimportant. One can then take over exact results available for the 
one dimensional Hubbard model in the dilute limit (see sec. |4.3.1|) . In particular. 



if the interactions added to ( |3.19| ) are repulsive (for Up = this corresponds to 
Kp < 1/m^), one has Up oc \n — 2/m|, but now K* = l/(2m^), and in particular 
Mo- oc (n — 2/m)^, i.e. there is a very strong divergence of the spin susceptibility. On 
the other hand, for attractive extra interactions in ( p.l9|) , a spin gap opens close to 



n = 2/m, and one then has Up oc \n — 2/m\, K* = \jvp? . A summary of the different 
types of critical behavior of the Luttinger liquid parameters Kp and Up^^j and of some 
derived physical quantities is given in table % 

The differences between even and odd filling fraction m may a priori seem 
surprising. However, there is a simple explanation, illustrated in figure |l]. Assume 
a lattice model with strong finite-range repulsion. Then, for m even, the ground 
state will essentially be a regular sequence of atoms, with one fermion every m/2 
lattice sites. The spins of the fermions will interact via an exchange term that is a 
generalization of the usual t^/f/ term, so that in fact we have a Akp charge density 
wave forming an antiferromagnetic spin chain. This is exactly the ground state we 
derived from ( |3.16| ). 



On the other hand, for odd m a completely equidistant arrangement is not 
possible, as illustrated in fig.|l] for m = 3. In fact, a band filling of 1/m corresponds 
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Table 1: Critical behavior of the different metal-insulator transitions considered in 
this paper. The transition is approached from the metallic side. For additional 
explanations on the last column see the text. 





n = 2/m 
m. even 


n = 2/m 
m. odd 


n 2/m 
Kp = const. 
m. even 


n 2/m 
Kp = const, 
m odd 


k; 




3/m2 


2/m2 


l/(2m2) 


Up 


const. 


const. 


\n — 1/m 


\n — l/m\ 


Ua 


const. 


const. 


const. 


{n — 


7 


const. 


const. 


1/ 72 — 1/m 


1/ n — 1/m 


X 


const. 


const. 


const. 


l/{n - l/mf 


Rw 


const. 


const. 
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to an average 2/m fermions per site. Consequently, the ground state arrangement is 
an alternation of short and long "bonds", and one has a 2kp CDW. In addition, there 
is of course an alternation of the exchange constants from long to short bonds, and 
such an alternation is well-known to introduce a gap in the spin excitation spectrum. 
We thus again explain the behavior derived from ( 3.17|) . 



The above considerations can be extended to obtain a lower bound on the 
correlation exponent Kp in a certain class of models: consider a lattice model with a 
finite range interaction of the form 

Hint = VmniUi+m ■ (3.20) 

i,m>0 

Let us assume further that Kn = for m > mo, where mo thus specifies the range 
of the interaction (mo = 1 and 2 for the Hubbard and extended Hubbard models 
discussed below). I will further assume strictly repulsive interactions, e.g. > Kn+i. 
Then, in the strong coupling limit it is quite clear that the longest interparticle 
distance in a stable structure is mo. In particular, the most dilute stable structure 
has just one particle every mo sites, i.e. it has period mo, and commensurate states 
with period larger than mo are not stable (stability here implies the existence of 
a gap in the charge excitation spectrum). In the continuum limit, a structure with 
period say mo + 1 can be stabilized by a term cos( v/8(mo + l)0p) , which has scaling 
dimension Xm^+i = 2(mo + l)'^Kp. This term will generate a stable commensurate 
structure if it is relevant, i.e. if Xmo+i ^ 2. On the other hand, for it not to create a 
stable commensurate structure (as expected), one needs Xmo+i > 2. Thus, for purely 
repulsive models with interaction of the form ( |3.20| ) and range mo one has the limit 



The assumption we have made here is that quite generally Kp decreases if the inter- 
action term ( p.20| ) is increased by applying a overall scale factor that is larger than 
unity. The limit ( |3.21| ) is satisfied for the lattice models discussed in the next section. 
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4 The Hubbard model in one dimension 



4.1 The Hamiltonian and its symmetries 

The Hubbard model is the prototypical model used for the description of correlated 
fermions in a large variety of circumstances, ranging from high-Tc superconductors to 
heavy fermion compounds and organic conductors. In spite of its apparent simplicity, 
there is still no general solution, or even a consensus on its fundamental properties. 



Notable exceptions are the cases of one and infinite dimensions 0]. In particular, 
in one dimension an exact solution is available. This solution gives exact energies of 
the ground state and all the excited states in terms of the solution of a system of cou- 
pled nonlinear equations. On the other hand, the corresponding wavefunctions have 
a form so complicated that the explicit calculation of matrix elements, correlation 
functions and other physical quantities has remained impossible so far. In the fol- 
lowing sections I shall describe in some detail the energy spectrum obtained from the 
exact solution. Subsequently, I will show how the knowledge of the energy spectrum 
can be combined with the results of the preceding two chapters to obtain a rather 
detailed picture of the low-energy properties, in particular of correlation functions, 
and of the metal insulator transition. 

The Hamiltonian in one dimension has the well-known form 

H = -tY,{clsCi+i,s + cl+i,sCi,s) + UY,ni,^ni^i , (4.1) 

i,s i 

where Cj_s is the fermion annihilation operator on site i with spin s, rii^g is the cor- 
responding number operator, and the sum is over the L sites of a one-dimensional 
chain with periodic boundary conditions. 



The model has two global SU(2) symmetries the first is the well- 

known spin rotation invariance, with generators 

C = E 4tc.a, = C. = \ j:{n.,i - n,t) . (4.2) 

i=i i=i 

The second type of symmetry is particular to the Hubbard model and relates sectors 
of different particle numbers. Its generators are 

V = E(-l)'ci>TCia' = ^^ = 2 ^('^ia + ^^.t) - ■ (4-3) 

1=1 i=l 

The total symmetry thus is SU{2) x SU{2) ^ 5*0(4). One should notice that more 
complicated interactions, e.g. involving further neighbors, will conserve the spin 
rotation invariance but in general not the "charge" SU(2) invariance ( |4.3| ). Rather, 
this second symmetry will become the standard global U{1) invariance associated 
with particle number conservation. It is nevertheless possible to construct particular 
types of further-neighbor interactions which do conserve the full SU{2) x SU{2) 
invariance. 
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4.2 The exact solution: ground state and excitations 

The exact wavefunctions of the one-dimensional Hubbard model for particles are 
superpositions of plane waves characterized by a set {kj} of momenta 0. The allowed 
values of kj are obtained from the solution of the coupled set of nonlinear equations 

^ ne( "^""7r"°M (4-4) 

(4.5) 





Here M is the number of down-spin electrons (M < N/2) and e{x) = {x + i)/{x — i). 
The Aq, are parameters characterizing the spin dynamics. We note that in general, 
both the fcj's and the A's are allowed to be complex. The energy and momentum of 
a state are 

N N 

E = -2tJ2 cos kj , P = J2 ■ (4-6) 
4.2.1 Solutions of the Bethe ansatz equations 

The determination of all the solutions of eqs. (|4.4|, [4.5|) is not easy. It has recently 



been shown (under certain assumptions) that these equations do indeed give all the 
"lowest weight" (with respect to SU{2) x SU{2)) eigenstates of the Hubbard model, 
i.e. all states satisfying ri\ip) = (\ip) = 0. The complete set of eigenstates then is 
obtained acting repeatedly with r^"'" or ("^ on {ip) Here I will limit myself to the 
ground state and to the low-lying elementary excitations. These questions have been 
investigated in some detail [|^, 0, however, the finite chain data presented below 
seem to be quite useful in understanding the nature of the excitations, and I therefore 
discuss them in detail. 

If both the fc's and the A's are all real, only the phases in ( [4.4| , [4.5| ) have to be 
determined. Taking the logarithm of these equations, one finds 

M 

L^j = 2TxIj + 2 ^ arctan[4(A« - sin kj)/U] (4.7) 

N M 

2^arctan[4(A„ - sin%)/f/] = 2TxJa + 2 ^ arctan[2(A„ - \p)/U] . (4.8) 

j=l /3=1 

The quantum numbers {/j} are all distinct from each other and are integers if M is 
even and half-odd integers (HOI, i.e. of the form 1/2, 3/2, . . . ) if M is odd, and 
are only defined modulo L. Similarly, the set {J a] are all distinct and are integers if 
— M is odd and HOI if — M is even. Moreover, there is the restriction 

|J„| < (A^-M + l)/2 . (4.9) 
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Summing ( [4.7| ) over j and ( ^.Sp over a, the total momentum is found as 



N M 



Ground state. The ground state is nondegenerate only if N is of the form 4i/ + 2 {ly 
integer): obviously, if N is odd, the ground state has (at least) spin 1/2. Further, if 
N is an integer multiple of 4 the noninteracting ground state has a sixfold degeneracy, 
and in the interacting case the ground state turns out to be a spin triplet. In the 
following I shall restrict myself to the case of N = Nq = Au + 2, i.e. the ground 
state is nondegenerate (in the following, A^^o will denote the particle number in the 
ground state). The ground state then is a singlet, with M = i.e. M is odd. 

The allowed values of the J's range from —{Nq/2 — l)/2 to {Nq/2 — l)/2. There are 
exactly Nq/2 such integers, i.e. all the J's are fixed. The J's are consecutive between 
-(A^o - l)/2 and {Nq - l)/2, i.e. 

{/,} = {-(iVo-l)/2,...,(iVo-l)/2} , (4.11) 
{Ja} = {-(iVo/2-l)/2,...,(iVo/2-l)/2} . (4.12) 

In the thermodynamic limit L ^ oo the distance between consecutive fc's or A's 
decrease like 1/L, and one can then find linear integral equations for the density 
of A;'s and A's on the real axis. Numerical results for the ground state energy as a 
function of particle density and U have been given by Shiba [59 . 



At half-filling, there is a gap in the charge excitations 0, and only spin ex- 
citations ( the ''2kp" triplet and singlet states below) are gapless. In particular, a 
finite energy is required to add or take out a particle. On the other hand, away from 
half-filling, both charge and spin excitations are massless. 



"4kF" singlet states. Excited states are obtained by varying the quantum num- 
bers. The first possibility, giving rise to excited singlet states, is obtained by removing 
one of the J's from the ground state sequence ( [4.11| ) and adding a "new" J: 

{/,} = {-(Aro-l)/2,...,-(Aro-l)/2 + zo-l, 

-(A^o - l)/2 + ^0 + 1, . . . , (A^o - l)/2, Jo} , (4.13) 
{Ja} = {-{No/2 -l)/2,..., {No/2 - 1) /2} , 

where |Jo| > (A'q — l)/2. This is a two-parameter family of excited states, called 
(somewhat misleadingly) "particle-hole excitation" by Coll. To understand the ex- 
cited states, we shall in the following consider systems of finite size, rather than taking 
the thermodynamic limit directly. It should however be quite evident how spectra 
like those of fig.^] develop into a true continuum in the thermodynamic limit. In fig.|^ 
we show numerical results for the energy-momentum spectrum of the states ( [4. 131 ) 
for a chain of 40 sites. The same states with k —k are obtained using negative 
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Iq. One notices a sharp minimum in the excitation energy at k/n = 1.1(1.7) = 4k p 
(this is why I call these excitations ^^kp" singlets). In the thermodynamic limit, the 
gap at Akp vanishes. These excitations are at the origin of the power-law behav- 
ior in the density-density correlation function (|2.47| ) around 4k p. Moreover, in the 



bosonization formalism, the Akp density correlations are entirely determined by the 
charge (0p) modes. Consequently, we identify the charge velocity Up as the slope of 
the excitation spectrum of fig.^ at = 0. Finally, we notice that the total number 
of states in this branch decreases as one approaches half-filling, and that this branch 
disappears altogether at half-filling. In fact, at half-filling, charge excitations have a 



gap and are described by solutions with complex /c's |52 



"2kF" triplet and singlet states. Excitations of the J's with all A's and A;'s real 
are only possible if M < N/2. The simplest excitations of this type are obtained 
considering M = N/2 — 1 which has total spin S = 1 (triplet). Now the restriction 
( [4.9|) allows N/2 + 1 different J's, i.e. we have two free parameters (the "holes" in 
the J-sequence), leading to sequences of quantum numbers of the form 

{/,} = {-iVo/2 + l,...,iVo/2} , 
Jl = -iVo/4 + 5„,,i , (4.14) 

Ja = Ja-1 + ^ + Sa,ai + Sa,a2 (a = 2, . . . , M) , 

where 1 < ai < ^2 < M, and Sa^p is the usual Kronecker symbol. Numerical results 
for this type of excitations are shown in fig.^ Corresponding states with negative k 
are obtained shifting the {J,} in ([4.14|) by one unit to the left. There now is a sharp 
minimum at k/n = 0.55 = 2kF, the gap again vanishing in the thermodynamic limit. 
In the long- wavelength limit, these states are the only spin-carrying excitations at 
constant particle number, so that the slope of the spectrum at A; = is equal to 
the spin velocity of the bosonized model. The low-energy excitations around 2kF 
are responsible for the spin contribution to the 2k p spin-spin correlations. As in 
fig.^ the structure of the excitations doesn't change much between weak and strong 
correlations, however, the energy scale does. In fact, the lowering of the energy scale 
in going to strong correlations corresponds to the lowering of the exchange energy 
(?s! At^ /U) in the strong correlation limit. The results of fig.|^ show apparent gaps at 
2kp and Akp. These are finite size effects: in the thermodynamic limit L —>■ oo the 
gaps vary like 1/L, and simultaneously the spectra develop into a continuum. 

Together with the triplet excitations ( [4. 141 ) there are also singlet states (M 



No/2), which are obtained by having one pair of complex conjugate A's among the 
solutions to the original equations ( [4.4|, [4.5|) . The energies of these states are shown 



by triangles in fig.^ It is remarkable that these states are nearly degenerate with the 
triplet states and in fact become exactly degenerate in the thermodynamic limit. 

The existence of singlets and triplets with the same energy shows that these 
states are in fact the combination of two noninteracting spin-1/2 objects, commonly 
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called spinous. Of course, because of total spin conservation, these objects can only 
be excited in pairs as long as one keeps the total number of particles fixed. 



Added particle. Adding one particle to the 4z/ + 2 ground state and leaving M 
unchanged, the J's and the J's are HOI. There are now M + 1 allowed values for the 
M distinct J's, and the low-energy states then are parameterized by 

{/,} = {-(Ao-l)/2,...,(Ao-l)/2,/o} , 
Ji = -M/2 + 5„,,i , (4.15) 

Ja = Ja-1 + I + Sa,ai (« = 2, . . . , M) , 

where |/o| > {Nq — l)/2, and 1 < ai < M. Corresponding spectra for different 
bandfillings are shown in fig.§ for Jq > 0. The symmetric spectra with negative k are 
again obtained using Jq < 0. The state of minimal excitation energy has momentum 
kp, as in the noninteracting case. The shallow arches in fig.§ correspond to varying 
ai, and are in fact of the same shape as the lowest branch (from k = 0.057r to 
k = 0.557r) in fig.|^, i.e. they correspond to single-spinon excitations. Close to kp 
the energy of theses states varies as u^{k — kp). On the other hand, varying Iq, one 
goes from one arc to the next, and the corresponding excitation energy (the upper 
limit of the quasi-continuum) varies as Up{k — kp). One also sees that going from one 
arc to the next in fig.^ (i.e. increasing Iq) the shape of an individual arc is basically 
unchanged. Varying Iq corresponds to a variation of the momentum of the added 
particle, and the figure thus shows that the total energy of a state is just the sum of 
the spinon energy and the "charge" energy associated with the added particle. One 
thus sees that charge and spin degrees of freedom do not interact. This is certainly in 
agreement with the predictions of the bosonization formalism, however, the fact that 
spin and charge separate even in highly excited states is special to the Hubbard model 
(a priori, the bosonized theory can only be expected to be an effective low-energy 
theory for the Hubbard and other lattice models). 

Another notable fact in fig.^ is the number of available states as Iq is varied: 
for Nq = 14, 22, and 30 there are respectively 13, 9, and 5 spinon arches. This means 
that without exciting the spins there are 26, 18, and 10 states for the extra particle 
available (counting states at negative k), i.e. the Iq branch stops at = vr — kp, 
rather than at A; = vr as in a noninteracting system. The explanation of this fact 
is rather straightforward for large U when double occupancy of sites is forbidden: 
in a system with L sites and Nq electrons, there are only L — Nq sites available at 
low energies. These states then form the "band" in fig.|[ There are of course states 
involving doubly occupied sites, however these are separated from the continuum of 
fig.^ by a gap (these states are solutions of ( [1.4| , O ) with complex /c's [52|). This 



separation of states occurs for any, even very small U, and can actually be shown by 
a perturbative argument |^ . 
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Added hole. Finally, let us consider states with one hole in the 4z/ + 2 groundstate 
(N = No - 1,M = {N - l)/2). Then both the /'s and the J's are integers. The 
energy is minimized chosing consecutive J's between — (A^ — l)/2 and (A^ — l)/2, but 
there are M + 1 possibilities for the M J's. States corresponding to the sequence 

{/,} = {-(iV-l)/2,...,(iV-l)/2} , 
Ji = -M/2 + 5,,,1 , 

Ja = Ja-l + l+Sa„a (a = 2, . . . , M) (4.16) 

are shown as the lowest arc between k = — 0.257r and k = 0.257r. This is a one-spinon 
branch, this state having necessarily S = 1/2, with velocity Mo-. One can of course 
also create a hole in the sequence of /'s. The energy spectrum for the quantum 
numbers 

h = -(iV + l)/2 + 5,,i 

= /,_i + l + (5,„, (j = 2,...,iV) (4.17) 

Ji = -M/2 + 5„,,i , 

Ja = Ja-l + l+6a„a (« = 2, . . . , M) (4.18) 

(where the free parameters obey 1 < ji < A^, 1 < ai < M) is also shown in fig.^. 
Similarly to fig.^, one notices that varying the "charge" quantum number ji one 
creates a branch with velocity Up. 



4.3 Low energy properties of the Hubbard model 
4.3.1 Luttinger liquid parameters 

In a weakly interacting system the coefficients Kp and Wjy can be determined pertur- 
batively. For example, for the Hubbard model one finds 

Kp = l-U/{iTVF) + ... , (4.19) 

where vp = 2tsin(7™/2) is the Fermi velocity for n particles per site. For larger U 
higher operators appear in the continuum Hamiltonian ( p. 301 ), e.g. higher derivatives 
of the fields or cosines of multiples of \/8(j)rj. These operators are irrelevant, i.e. they 
renormalized to zero and do not qualitatively change the long-distance properties, 
but they do lead to nontrivial corrections to the coefficients Ui,, Kp. In principle 
these corrections can be treated order by order in perturbation theory. However, this 
approach is obviously unpractical for large U, and moreover it is at least possible that 
perturbation theory is not convergent. To obtain the physical properties for arbitrary 
U a different approach is therefore necessary. 

I note two points: (i) in the small-?/ perturbative regime, interactions renor- 
malize to the weak-coupling fixed point gl = 0, K* = 1; (ii) the exact solution 
does not show any singular behavior at nonzero U, i.e. large U and small U are the 
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same phase of the model, so that the long-range behavior even of the large U case is 
determined by the fixed point gl = 0. Thus, the low energy properties of the model 
are still determined by the three parameters Up^^ and Kp. 

The velocities be obtained from the long wavelength limit of the '%kp^^ 

and excitations discussed above. In the thermodynamic limit the corresponding 

excitation energies are easily found from the numerical solution of a linear integral 



equation jSQ]. Results are shown in fig.y for various values of U/t. Note that for 
U = one has Up = = 2t sin('7rn/2), whereas for U ^ oo Up = 2tsin(7rn), 
Ua = (27rt^/f/)(l — sin(27m)/(27rn)). In the noninteracting case oc n for small 
n, but for any positive U oc in?. The Wilson ratio, eq.( |2.38| ), obtained from the 
velocities is shown in fig.^ For [7 = one has Rw = 1, whereas for ?7 — > oo Rw = 2 
for n 7^ 1. 

To obtain the parameter Kp from the exact solution note that the gradient of 
the phase field 0p is proportional to the particle density, and in particular a constant 
slope of (j)p represents a change of total particle number. Consequently, the coefficient 



Up/Kp in eq. (|2.31|) is proportional to the variation of the ground state energy Eq 



with particle number IT 

1 d^E^{n) _'K Up _ 1 ^^20) 



L din? 2 Kp Tn?K 



Equation ( [4.20| ) now allows the direct determination of Kp-. Eo[n) can be obtained 



solving (numerically) Lieb and Wu's integral equation, and Up is already known. 
The results for Kp as a function of particle density are shown in figj^ for different 
values of U/t. For small U one finds in all cases agreement with the perturbative 
expression, eq. ( ^4.19 ), whereas for large U Kp 1/2. The limiting behavior for 



large U can be understood noting that for U = oo the charge dynamics of the system 
can be described by noninteracting spinless fermions (the hard-core constraint then 
is satisfied by the Pauli principle) with kp replaced by 2/?^?. Consequently one finds a 
contribution proportional to cos{4kpx)x~'^ in the density-density correlation function, 
which from eq. ( 2.47|) implies Kp = 1/2. One then finds an asymptotic decay 



like cos(2A;^x)x In ^^(x) for the spin-spin correlations, eq.( p.48|) , and an exponent 



5 = 1/8 in the momentum distribution function. The result 6 = 1/8 has also been 
found by Anderson and Ren |5^, and by Parola and Sorella |5^. Ogata and Shiba's 
numerical results are quite close to these exact values. 



As is apparent from fig.^, the strong-coupling value Kp = 1/2 is also reached in 
the limits n — 0, 1 for any positive U. For n ^ this behavior is easily understood: 
the effective interaction parameter is U/vp, but vp goes to zero in the low-density 
limit (corresponding to the diverging density of states). The limit n — >■ 1 is more 
subtle: in this case nearly every site is singly occupied, with a very low density of 
holes. The only important interaction then is the short range repulsion between holes, 
which can be approximated by treating the holes as a gas of spinless noninteracting 
fermions. Using (|4.2CI| ), one then again finds Kp = 1/2. 
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We note that in the whole parameter region, as long as the interaction is repul- 
sive one always has Kp < 1, which means that magnetic fluctuations are enhanced 
over the noninteracting case. On the other hand, superconducting pairing is always 
suppressed. 

The results of fig.^] are valid for n — > 1, but not for n = 1. In the latter case, 
there is a gap in the charge excitation spectrum, as expected from the umklapp term 
( |3.2|) , and the correlations of 0p become long ranged {Kp^^ff = 0). Close to half-filling, 
the asymptotic behavior of correlation functions then exhibits a crossover behavior 
from half-filled-like behavior at short distances to the general form at long distances, 
as discusses in secIO] above (see in particular eq. ( |3.10| ). 



Results equivalent to the present ones can be obtained using the conformal 
invariance of the Hubbard model 0. These results have subsequently be gener- 



alized to the case with an applied magnetic field |58 



4.3.2 Transport properties 

The exact solution of Lieb and Wu can also be combined with the long-wavelength 
effective Hamiltonian ( p.3(]| ) to obtain some information on the frequency-dependent 



conductivity a{u!). On the one hand, from eq. (|2.39| ) there is a delta function peak 
at zero frequency of weight 2KpUp. On the other hand, the total oscillator strength 



is proportional to the kinetic energy |59]: 



<^tot 



/oo 
a{tu)du = -n{Hkin)/L . (4.21) 
-oo 



Thus, both the weight of the dc peak and the relative weight of the dc peak in 
the total conductivity can be obtained and are plotted in fig.|^. As expected, far 
from half-filling, all the weight in atot is in the dc peak. For exactly half-filling the 
dc conductivity vanishes, due to the existence of a gap Ac for charge excitations 
created by umklapp scattering, and all the weight is at > Ac. Fig. 2 then shows 
that as n ^ 1 umklapp scattering progressively transfers weight from zero to high 
frequency. The crossover is very sharp for small or large U, but rather smooth in 
intermediate cases {U/t ^ 16). This nonmonotonic behavior as a function of U can 
be understood noting that initially with increasing U umklapp scattering plays an 
increasingly important role. Beyond U /t ^ 16, however, the spinless-fermion picture 
becomes more and more appropriate, and sX U = oo one again has all the weight 
in the dc peak. The linear vanishing of cxo as n — 1 implies a linear variation of 
the ratio n/m* with "doping". By the thermopower-argument of sec. p.l| carriers 
are hole-like for < 1, and electron-like for n > 1, provided one is close to the 
transition. The thermopower of the one-dimensional Hubbard model as well as the 
rather subtle crossover occuring in the vicinity of the critical point n = 1, [/ = have 
been analyzed in detail recently |6T1 . 
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4.3.3 Spin charge separation 

The Hubbard model also provides a rather straightforward interpretation of the spin- 
charge separation discussed above. Consider a piece of a Hubbard chain with a half- 
filled band. Then for strong U there will be no doubly-occupied sites, and because 
of the strong short-range antiferromagnetic order the typical local configuration will 
be 

■■■nnnnnn--- 

Introducing a hole will lead to 

■■■nni onnn--- 

and after moving the hole one has (note that the kinetic term in the Hamiltonian 
does not flip spins) 

■■■nonninn--- 

Now the hole is surrounded by one up and one down spin, whereas somewhere else 
there are two adjacent up spins. Finally, a few exchange spin processes lead to 

■■■nonnnin--- 

Note that the original configuration, a hole surrounded by two up spins has split into 
a hole surrounded by antiferromagnetically aligned spins ("holon") and a domain- 
wall like configuration, two adjacent up spins, which contain an excess spin 1/2 with 
respect to the initial antiferromagnet ("spinon"). The exact solution by Lieb and Wu 
contains two types of quantum numbers which can be associated with the dynamics of 
the spinous and holons, respectively. We thus notice that spinous and holons [^, ^ 
have a well-defined meaning in the present one-dimensional case. 

The above pictures suggest that, as far as charge motion is concerned, the Hub- 
bard model away from half-filling can be considered as a one-dimensional harmonic 
solid, the motion of the holes providing for an effective elastic coupling between ad- 
jacent electrons. This picture has been shown to lead to the correct long-distance 



correlation functions for spinless fermions |6J, For the case with spin, this sug- 



gests that one can consider the system as a harmonic solid with a spin at each site 
of the elastic lattice (lattice site = electron in this picture). 

Let us now show that this gives indeed the correct spin correlation functions. 
In a continuum approximation, the spin density then becomes 

cTix) = J2SmSix-x^) , (4.22) 

m 

where the sum is over all electrons. After a Fourier transformation of the delta 
function the spin-spin correlation function becomes 

(cTix) ■ o-(O)) = -i-^ J dqdq' J2 e-'^^Sra ■ 5„,,e^('?^™+'?'^™')) . (4.23) 

' m.m' 
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The exchange energy between adjacent spins is always antiferromagnetic, whether 
there is a hole between them or not, and consequently the low-energy spin dynamics 
always is that of an antiferromagnetic chain of localized spins. Under the additional 
assumption that the spin-spin correlations on an elastic lattice depend mainly on the 
average exchange constant and not so much on the fluctuations induced by motion of 
the electrons, the average in ([4.23| ) factorizes into separate spin and charge factors. 



Following the hypothesis about harmonic motion of the electrons, we write Xm = 
Rm+Um, where Rm = m/\l—n\ is the average position of the mth electron and Um the 
displacement with respect to this position. Note that the "harmonic solid hypothesis" 
implies that Um+i — Um is small, but not necessarily Um and Um+i separately. In the 
averages over atomic positions now all terms with q ^ q' vanish, and one has 



cr X 



) ■ o-(O)) ^ JdqJ2 e-^'^^Sm ■ 5^.)e*^(^™-^-')(e*'("'"-"-')) . (4.24) 



m,m' 



The average over Um in ( [4.24|) has a power law behavior: 



e 



with a{q) oc g^, i.e. it has a smooth g-dependence. On the other hand, the long- 
distance behavior of the spin-spin correlations of an antiferromagnetic spin chain is 
63, |33|, 131 



[Sm ■ S^,) ^ (-l)"'|m - m'\~' ln^/2 



m — m 



Therefore in ([4.24|) the g-integration is dominated by terms with g f» 7r(l — n) = 2ki 
Replacing the weakly g-dependent exponent a(g) by a{2kF) one obtains 



o-(O)) ^ y" c/ge-^"^ ^ (5„ ■ 5„/)e*'?(^™"^-' V - ^T^^*"^^ (4-25) 

m,m' 

= cos(2fci.x)x-^-"(2'=^Mn^/2^x) . (4.26) 



With the identification a{2kF) = Kp, this is precisely the result ( 2.48| ). We thus have 



shown that the spin-spin correlations of a correlated electron system can in fact be 
understood as those of an elastic lattice of spins. In that picture, the motion of the 
holes then only provides the effective elasticity for the lattice. 



4.3.4 The metal— insulator transition 

The one-dimensional Hubbard model is insulating for n = 1, U > 0, but conducting in 
all other cases. Moreover, in agreement with the discussion of section ^ , the metal- 
insulator transition occurs in a different way according to whether the interaction 
strength is varied at constant carrier density or whether one varies the carrier density 
at constant U. 

It seems worthwhile here to compare the metal-insulator transition occuring 
as n — s> 1 with other scenarios for strongly correlated fermion systems in higher 
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dimension (see the review by Vollhardt [0). In the "nearly locahzed" picture, ef- 
fective mass effects predominate and enhance both the specific heat and the spin 
susceptibihty. Consequently, the Wilson ratio (1/(1 + Fq) in Fermi liquid language) 
remains nonzero as the metal-insulator is approached. On the other hand, in the 
"nearly ferromagnetic" (or paramagnon) picture, only the spin susceptibility is en- 
hanced significantly, and therefore Rw can be much larger than unity. The behavior 
found here in the one-dimensional case is quite different from both these scenarios: 
generally Rw < 2, and approaching the metal-insulator transition Rw — ^ 0. This 
occurs because generally an enhancement of the mass of the charge carriers (i.e. a 
decrease of Up) has no influence on the spin degrees of freedom (see fig.l). This is 
rather straightforwardly understood in terms of spin-charge decoupling, as explained 
in the previous section: charge and spin excitations move nearly independently of 
each other, and in particular the spin dynamics is determined by antiferromagnetic 
nearest-neighbor exchange. In particular the spin susceptibility remains finite even 
when the mass of the charge carrier approaches infinity. 

Let us discuss the metal-insulator transition in more detail. The fact that Up 
and (Jo vanish linearly as n — > 1 seems to be consistent with a divergent effective 
mass at constant carrier density because Up ~ l/m*,ao ~ n/m*. A constant carrier 
density is also consistent with the fact that kp = 7m/2 is independent of U. It is 
not consistent with the hole-like sign of the thermopower as n ^ 1 from below, nor 
with the electron-like sign as n — >^ 1 from above: if the carriers are holes, the carrier 
density is the density of holes: n* = 1 — n. Treating the holes as spinless fermions, 
as already mentioned before, one expects (Tq — ^ because n* — 0, and 7 — > cxd 
because the density of states of a one-dimensional band diverges at the band edges. 



This agrees with what was found explicitly in section |4.3.1| . What is not so easily 



understood in this picture is the fact that kp (i.e. the location of the singularity of 
rifc) is given by its free-electron value 7rn/2, rather then being proportional to n*. 
One should however notice that is given by the single-particle Green's function 
which contains both charge and spin degrees of freedom. The location of kp then 
may possibly be explained by phase shifts due to holon-spinon interaction. This is 
in fact suggested by the structure of the wavefunction of the exact solution |jl7|. 

The magnetic properties do not agree with what one expects from an effective 
mass diverging as n — 1: and therefor x remain finite. Moreover, the NMR 
relaxation rate would have the behavior 1/Ti = aT+pVT, where the first (Korringa) 
term comes from fluctuations with q ^ 0, whereas the second term comes from 
antiferromagnetic fluctuations with q ^ 2kp. None of these properties is strongly 
influenced by the diverging effective mass observed e.g. in the specific heat. This 
fact is of course a manifestation of the separation between spin and charge degrees 
of freedom. 
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4.3.5 Other models 



It is clearly interesting to go beyond the Hubbard model. A simple generalization is 
the "extended Hubbard model" which includes a nearest-neighbor repulsion: 

= -t^(a|^ai+i,, + a|_^.i,ai^) + f/^raj-frai^ + \/^nir2i+i , (4.27) 

i,s i i 

For this model, exact eigenvalues can not be obtained in the thermodynamic limit. 
The parameters in eq. ( |4.20|) can however be calculated reliably for finite systems \dE\. 



In particular, at quarter filling one finds a Luttinger-liquid ground state, with Kp = 
1/4 at the metal-insulator transition which occurs with increasing V, in agreement 
with the discussion of sec. 

Exact exponents can be obtained for the model ([4.27|) in the limit f/ — > oo: then 
one has effectively spinless fermions (with kp 2kp) with nearest neighbor interac- 
tion, a model which can be exactly solved using the Jordan- Wigner transformation 
into the XXZ spin chain. In particular, the 4A;j7-component of ( p. 47] ) is related to the 



correlation function of Sz- From the known results one obtains, for a quarter 



filled band {n = 1/2), Kp = 1/(2 + (4/7r) sin-^(t;)), Up = irtVl^v^ / cos~\v) , with 
V = V/2\t\. Now Kp < 1/2 is possible. For v > 1 the system is in a dimerized insulat- 
ing state. Approaching the insulating state from f < 1 both Kp and Up remain finite, 
i.e. o"o jumps to zero at f = 1. For n ^ 1/2 the parameters Up,Kp can be obtained 
from numerical results ||39|. Quite generally, one has Kp > 1/8, but Kp = 1/2 for 



n — *• 0, 1, independent of v. On the other hand, — ^ as n — > 1/2 for t> > 1, i.e. 
in that case the weight of the dc conductivity goes to zero continuously, the point 
(f,n) = (1, 1/2) is thus highly singular. The same type of singularity also occurs at 



U = 0,n = 1 in the Hubbard model |60l]. Interestingly enough, one has Kp > 1 if 



V < -V2\t\, i.e. a finite amount of nearest-neighbor attraction is sufficient to lead 
to divergent superconducting fiuctuations even for infinite on-site repulsion. Also 
note that the singularities in Up and Kp at f = — 1 (attractive interaction) represent 
a point of phase separation. 

For the t — J model, there is one exactly solvable point (t = J) where exact 
exponents can be found using the Bethe ansatz [^. Away from this point, eq. 



( |4.20D has been used to obtain Kp from numerical data |70l. It is not easy to study 



the metal-insulator transition occuring for n — > 1 numerically, however the results are 
consistent with Kp = 1/2 in this limit. For large J there is a phase with predominantly 
superconducting fluctuations {Kp > 1). 



5 Conclusion 

In this paper, we have seen that using the bosonization method and exact solutions, 
one obtains a rather complete picture of many physical properties of interacting 
one-dimensional fermions. Probably the most important feature arising is the Lut- 
tinger liquid like behavior, characterized by non-universal power laws, together with 
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the separation of the charge and spin dynamics. One should also notice that there 
are no qualitative differences between weak and strong correlation. Metal-insulator 
transitions occur for commensurate bandfillings, and in particular at half-filing, for 
repulsive interactions. The transitions at varying particle density show qualitatively 
different behavior according to whether the bandfilling is an even or odd fraction. 

Another type of metal-insulator transition, not discussed in detail here, oc- 
curs in the presence of disorder (see ref. [0 and references therein). In this case, 
as is well-known, in the absence of interactions all states are localized. Repulsive 
interactions enhance localization. On the other hand, sufficiently strong attraction 
can lead to delocalization. In cases where there is no spin gap, the transition to 
the delocalized conducting (in fact superconducting) state occurs at Kp > 2, and at 
the transition one has a non-universal value of the correlation exponent, satisfying 
only K* > 2. In the case with spin gap, the superconducting state only occurs for 
Kp > 3. In this case, at the transition one has a universal value K* = 3. One should 
notice that these large values of Kp in fact correspond to rather strong attractive 
electron-electron interactions, and it seems doubtful that this can be realized in an 
experimental system. 

Unambigous experimental observation of Luttinger liquid like behaviour and of 
the associated metal-insulator transitions is made difficult by the fact that all possi- 
ble known candidates are in fact only ^itasz-one-dimensional, e.g. they consist of a 
parallel arrangement of conducting chains. At sufficiently low temperatures one thus 
always crosses over into a regime of effectively three-dimensional behaviour, charac- 
terized in particular by the occurence of different kinds of ordered states. The specif- 
ically one-dimensional behaviour is thus not always easily identified. Nevertheless, a 
number of very interesting experiments do exist. An early and spectacular example 
is the observation of diffuse X-ray scattering at wavevector 4k p in the compound 
TTF-TCNQ In fact, in a perturbative theory, no such scattering is expected. 



However, within the boson picture of the Luttinger liquid, such scattering is indeed 
expected for strongly repulsive interaction |j72|, ^ (see also eq. ( p.47| ) and the subse- 



quent discussion). More recently, NMR data on the series of compounds (TMTSF)2X 



have shown Luttinger liquid like behavior [74|, and in particular power-law depen- 



dence of the relaxation rate on temperature. In some of these compounds umklapp 
scattering is sufficiently strong to induce a (relatively small) Mott-Hubbard insulat- 
ing gap, and one can then observe a crossover between high-temperature metallic 
and low-temperature insulating behavior at a temperature between lOOK and 200K. 
The most recent observation concerns photoemission on (TMTSF)2PF6 0, where the 
spectral density does not show a Fermi edge like behavior, but is rather reminiscent 
of a power-law behaviors as expected in one dimension (see eq. ( |2.42|) ). However, the 
observed exponent is 5 ~ 1, a value that would imply very strong electron-electron 
repulsion. One then would expect strong effects of umklapp scattering, i.e. typically 
a non-metallic conductivity. However, experiment shows good metallic behavior in 
the temperature region concerned. The photoemission results thus do not seem to be 
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fully understood. It should be pointed out here that the metal-insulator transitions 
observed always occur at constant band-filling, upon varying temperature, pressure, 
or chemical composition (which in the last two cases is equivalent to changing the in- 
teraction strength). Unfortunately, it has up to now been impossible to dope organic 
conductors in a sufficiently controlled way that would make the transition occuring 
as a function of doping observable. 

Anderson has suggested that Luttinger liquid behavior might also occur in two 
dimensions |^5| , as well as in coupled chain systems . Under which circumstances 



this suggestion is correct does not currently clear. At least for the coupled-chain case, 
Anderson's suggestion is in contradiction with standard scaling |]7^ and renormal- 



ization group arguments which indicate that interchain hopping is a strongly 

relevant perturbation and therefore most likelily will destroy the Luttinger liquid 
behavior. 
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Figure 1: Charge arrangement for a system with strong on-site and nearest-neighbor 
repulsion for a third-filled (top) and quarter-filled (bottom) band. Large and small 
dots are occupied and empty sites, respectively. Note the alternation of nearest- 
neighbor distances in the third-filled case. 



Figure 2: "Akp" singlet excitation spectrum for a Hubbard chain of 40 sites with 22 
and 34 electrons. The lowest "arc" from k/n = 0.05 to k/n = 1.1 (or k/n = 1.7) 
is obtained by varying at fixed Iq = (Nq + l)/2 (cf. eq. the higher arches 

correspond to increasing Iq up to [L — l)/2. 



Figure 3: '"Ikp" spin singlet (A) and triplet (o) excitation spectrum for Hubbard 
chains of 40 sites with 22 electrons for different interaction strengths. 



Figure 4: Excitation spectra for one particle added into a Hubbard chain of 40 sites 
with 14, 22, and 30 electrons. The shallow arches correspond to varying ai (cf. 
eq. I \4-1^ ) at constant Iq, and Iq increases from one arc to the next. Zero energy 
corresponds to the A'^o + 1 particle ground state. 



Figure 5: Excitation spectra for an added hole in a Hubbard chain of 4O sites with 2^ 
electrons. The shallow arches correspond to varying ai (cf. eq. ( j.lS ) at constant ji, 
and ji increases from one arc to the next. 



Figure 6: The charge and spin velocities Up (full line) and (dashed line) for the 
Hubbard model, as a function of the band filling for different values of U/t: for 
U/t = 1, 2, 4, 8, 16 from top to bottom, for Up U /t = 16, 8, 4, 2, 1 from top to bottom 
in the left part of the figure. 



Figure 7: The Wilson ratio Rw for the one- dimensional Hubbard model, as a function 
of the band filling for different values of U/t (U/t = 16, 8, 4, 2, 1 for the top to bottom 
curves). 
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Figure 8: The correlation exponent Kp as a function of the bandfilling n for different 
values ofU (U/t = 1,2, 4, 8, 16 for the top to bottom curves). Note the rapid variation 
near n—1 for small U . 



Figure 9: Top: The weight of the dc peak in cr(a;) as a function of bandfilling for 
different values of U/t (U/t = 1,2, 4,8, 16 for the top to bottom curves). 
Bottom: Variation of the relative weight of the dc peak in the total conductivity 
oscillator strength as a function of the bandfilling n for different values ofU: U/t = 1 
(full line), 4 (dashed), 16 (dash-dotted), 64 (dotted), and256 (dash-double- dotted). 
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Figure 1: Charge arrangement for a system with strong on-site and nearest-neighbor 
repulsion for a third-filled (top) and quarter-filled (bottom) band. Large and small 
dots are occupied and empty sites, respectively. Note the alternation of nearest- 
neighbor distances in the third-filled case. 
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L = 40, N = 22, U/t = 16 
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Figure 2: "Akp" singlet excitation spectrum for a Hubbard chain of 40 sites with 22 
and 34 electrons. The lowest "arc" from k/n = 0.05 to k/n = 1.1 (or k/iv = 1.7) 
is obtained by varying io at fixed Iq = (Nq + l)/2 (cf. eq. (4-13)), the higher arches 
correspond to increasing Iq up to (L — l)/2. 
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Figure 3: "2kF" spin smglet (A) and triplet (o) excitation spectrum for Hubbard 
chains of 40 sites with 22 electrons for different interaction strengths. 
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L = 40, N = 14 +1, U/t = 16 
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Figure 4: Excitation spectra for one particle added into a Hubbard chain of 40 sites 
with 14, 22, and 30 electrons. The shallow arches correspond to varying ai (cf. 
eq. (4-15)) at constant Iq, and Iq increases from one arc to the next. Zero energy 
corresponds to the Nq-\-1 particle ground state. 
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L = 40,N = 22-l,U/t= 16 
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Figure 5: Excitation spectra for an added hole in a Hubbard chain of 40 sites with 22 
electrons. The shallow arches correspond to varying ai (cf. eq. (4-18) at constant ji, 
and ji increases from one arc to the next. 
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Figure 6: The charge and spin velocities Up (full line) and (dashed line) for the 
Hubbard model, as a function of the band filling for different values of U/t: for u„ 
U /t = 1, 2, 4, 8, 16 from top to bottom, for Up U /t = 16, 8, 4, 2, 1 from top to bottom 
in the left part of the figure. 
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Figure 7: The Wilson ratio Rw for the one-dimensional Hubbard model, as a function 
of the band filling for different values of U/t = 16, 8, 4, 2, 1 for the top to bottom 
curves). 
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Figure 8: The correlation exponent Kp as a function of the bandfillmg n for different 
values ofU (U/t = 1, 2, 4, 8, 16 for the top to bottom curves). Note the rapid variation 
near n = 1 for small U. 
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Figure 9: Top: The weight of the dc peak in cr(uj) as a function of bandfiUing for 
different values of U/t (U/t = 1,2,4,8,16 for the top to bottom curves). 
Bottom: Variation of the relative weight of the dc peak in the total conductivity 
oscillator strength as a function of the bandfiUing n for different values ofU: U/t = 1 
(full line), 4 (dashed), 16 (dash-dotted), 64 (dotted), and 256 (dash-double-dotted). 
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